Variational study of a two-level system coupled to a harmonic oscillator in a ultra-strong coupling 

regime 
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The non-classical behaviors of a two-level system coupled to a harmonic oscillator are investigated in the 
ultra-strong coupling regime. We revisit the variational solution of the ground state and find that the existing 
solutions do not account accurately for non-classical effects such as squeezing. We suggest a new trial wave 
function and demonstrate that it has an excellent accuracy on the quantum correlation effects as well as on 
energy. 
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A two-level system interacting with a harmonic oscillator 
appears in various fields in physics, ranging from an atom cou- 
pled to a photon inside of an optical cavity | ll to a cooper-pair 
box coupled to a nanomechanical oscillator ||2J. A theoretical 
model that describes the system is expressed by the Hamilto- 



man. 
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A (a + a') fTj, + 



(1) 



where ujq is the energy of the oscillator, is the spin level 
splitting, and A is the coupling strength. The Pauli matrices 
characterize the two level system, while a and denote the 
boson operators. 

The Hamiltonian ([T]) reveals completely different physics 
at different scales of wq, and A. For example, a cavity 
quantum electrodynamics (QED) operates in a regime where 
\,\ujQ — Sl\ <S£^ bjQ^^, under which the counter-rotating terms 
of the Hamiltonian ([I}, aV+ and acr^, can be neglected. 
Within this rotating wave approximation (RWA), the model 
reduces to a Jaynes-Cummings model and is solvable ex- 
actly |3|. The ground state is a simple direct product of the 
ground states of the oscillator and the spin. On the other 
hand, in the so-called ultra-strong coupling regime (il ^ 
and A ^ -v/Owo), there arise interesting quantum effects in 
the ground state; the two-level system and the oscillator are 
entangled. The degree of entanglement increases monoton- 
ically as a function of the dimensionless coupling constant 
g = 2A/v^ncjo. It also shows a squeezing effect m, that 
is, the variance of momentum quadrature in the ground state 
becomes smaller than the uncertainty minimum L5^-^. Here- 
after, we will set ujq — 1 for simplicity. 

The aforementioned ultra-strong regime has recently at- 
tracted much attention because of the possibility to realize 
it experimentally in a circuit QED system fT UTOl . Further- 
more, as rather a theoretical problem, the ultra-strong cou- 
pling regime when ft becomes oo has been also studied ex- 
tensively because it shows a phase transition-like behavior in 
the entanglement between the oscillator and the qubit in the 
ground state lfTTHT6]| . 



In this paper, we develop a variational ground state. The 
variational method had already been used for the Hamilto- 
nian ([TJ. A displaced coherent state for g ^ 1, a displaced 
squeezed state for g <C 1 ifTTl . and a superposition of two 
coherent states for 5 ~ 1 ifTSl [T9l have been suggested as 
a variational ground state, respectively. The accuracy of the 
suggested variational states were reasonable on the ground- 
state energy lfT7ljT9i . However, we find that the variational 
states suggested in the previous works |fT8l[T9l underestimates 
the squeezing effect significantly in the intermediate regime 
{g ~ 1); see Figs. [2] (b) and (e). The failure gets more pro- 
nounced as il increases. 

This observation implies that an energy-optimized vari- 
ational wave function does not necessarily capture all the 
quantum correlation effects in the true ground state of the sys- 
tem. In this paper, we suggest a new variational ground state 
which captures accurately the squeezing effect. 

First of all, we derive an equivalent model with only a boson 
degrees of freedom, followed by a derivation of an effective 
classical Hamiltonian by utilizing a parity symmetry in the 
model |20|. We observe that this effective classical Hamilto- 
nian shows bifurcation around 5 = 1, which is consistent with 
the observation in Ref. ( lfT2l[l4ll ) based on a classical model of 
the Eq. ([T]l, and gives us qualitative account for the properties 
of the ground state. We then suggest a superposition of two 
displaced-squeezed state as a trial wave function, and demon- 
strate that it can predict the squeezing effect accurately. This 
shows that a deformation of each superposed wave packet is 
crucial to the squeezing effect of ground state. We also dis- 
cuss the drastic change in the entanglement degree as g varies 
across g = 1 when il — > 00. 



I. EFFECTIVE HAMILTONIAN 

The model (fTl) has a useful symmetry. If \n' ^a') = 
H |n, a), then n -\- a' (mod 2) = n + a (mod 2), where 
|n, (t) = \n) (E) |cr), |n) (n = 0, 1, 2, • • • ) is the boson num- 
ber state, a'aln) — n\n), and \a) (cr = ±1) the spin state, 
CT^ |cr) — a\a-). We further stress this symmetry by introduc- 
ing a generalized "parity" operator 
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FIG. 1: (color online) The profile of H~^{x,0) (a) for Q = 10 and 
(b) for — 10^. Different curves correspond lo g = 0.1, 1,2,3 
(from bottom to top). The inset in (b) zooms in the region around 
x = 0. 



Clearly, the Hamiltonian ([T]) is invariant under transformation 
described by 11, [H, H] = 0. The operator H has two eigen- 
values ±1, and corresponding eigenstates are given by 



(at) 



\0,a) (a = ±l) 



(3) 



Namely, n \ (p^) = ± Ifn)- The Hilbert space can be decom- 
posed into a direct sum £ = of two subspaces £± 
spanned by respectively. Accordingly, the Hamiltonian 
can also be written as H — where H± belongs to 
the subspace £±, and the partition function Z = Ti[e^^^] as 
Z ^Z+ + where Z± = Tr[e-'3^±]. 

Now we make a key observation, suggested in the expres- 
sion for the parity eigenstates in Eq. (|3]l, that the combina- 
tion a)ax behaves like a boson operator within each sub- 
space £±. To show this rigorously, we first take a drone- 
fermion representation 121] of spin = 2/^/ — 1 and 
a+ = (tI_ = {d + d^), where / and d are fermion operators. 
We then define a new boson operator = [p — f){d + d^) 
which satisfies [b, 6^] = 1. The parity basis states are then 
written as 



(6t)" 



|0,a) 



(4) 



where |0, a) serves as the "vacuum" of the new boson operator 
b within £„. It is now clear that the sub-block Hamiltonians 
can be rewritten solely in terms of the boson operator b as 



(5t-A)(&- 



A) ± ^r2cos(7r6t6) , 



(5) 



The same expression was derived also in Ref. 1 201 using a 
canonical transformation. Having only a boson field, the ex- 
pression (|5]l allows us to compare the model directly with cor- 
responding classical system by deriving the classical effective 
potential. 

To get a classical effective potential, we express the par- 
tition functions Z± in terms of a functional integration ll22l . 
Z± = J 2?[0* (r), 0(t)] exp(—5'±) with the action defined by 



S+ - 



/3 



dT 



+ (0* - A)((/) - A) ± Jf7e~20> 



(6) 

where f3 is the inverse temperature, d-r denotes the partial 
derivative with respect to the imaginary-time r. Making a 




3 4 0.0 0.5 1.0 1.5 2.0 



FIG. 2: (color online) Comparison of variational solutions with the 
exact solution. The (blue) empty square represents the DCS wave 
function, while the (red) filled circle indicates the DSS wave func- 
tion. (a,e) Error in the energy with respect to the exact energy. (b,f) 
Error in the tangle. (c,g) Error in the momentum variance. (d,h) The 
momentum variance of the ground state, the (black) line indicates the 
exact solution. 



change of variables, p = (0* + 0)/a/2 and p ~ — 
(f))/V2, we can rewrite the functional integrals as Z± = 
JV[x{T),p{T)]exp{-S±) with 

S±= dT[-ipdrX + H^,{x,p)]. (7) 

Jo 

Here H^(x,p) is the effective classical Hamiltonian defined 
by 
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(8) 



The potential profile is depicted in Fig.[T]for both ft — 10 
and = 10^. For 5 <C 1, the classical potential has one 
local minimum and is highly anharmonic illustrating why a 
squeezed state localized around x = is a good variational 
solution in this regime. As g approaches the critical value, 
1, it shows a bifurcation of its local minimum as in the cor- 
responding classical system. As g increases further, its local 
minimum around the origin disappears and a single local min- 
imum develops at a; ^ gy/fl/2. It is also consistent with the 
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fact that the coherent state is a good approximation in this 
regime. 



II. VARIATIONAL SOLUTIONS 

Above we have seen from the effective classical potential 
that the model in Eq. ([T]l shows highly non-trivial behaviors 
in the regime g ^ 1. Now we investigate more closely this 
regime in terms of variational wave functions. To test the ac- 
curacy of a given trial wave function, we will examine not 
only the energy but also the quantum correlation effects such 
as squeezing and spin-oscillator entanglement. We character- 
ize the squeezing effect with the momentum variance AP and 
the spin-oscillator entanglement with tangle [23 1 defined by 
T — 2[1 — Tr(pg)], where ps is a reduced density matrix of 
the two-level system. We test the variational solution for two 
cases n = 10 and SI = 10^ by comparing it with the exact 
solution. The exact solution is numerically calculated by the 
exact diagonalization method ll24l . where a sufficient number 
of Pock basis states |n) (typically 300 states for SI = 10 and 
3000 states for SI — 10^) are kept until the desired accuracy is 
achieved. 

In Refs. ifTTlfTOl . a superposition of two displaced coherent 
states was suggested as a trial wave function, 



\^p{ai,a2,t)) = (1 -t)\ai) +t\a2), 



(9) 



where |a) — D{a) |0) with _D(q:) ~ exp[a(&^— 6)]. Hereafter 
we will refer to the trial wave function in Eq. (j9]l as the double 
coherent state (DCS). It shows a reasonable accuracy for the 
energy (with error about 1 %) in the intermediate regime [see 
Figs. [2] (a) and (e)]. Although neither our effective classical 
potential nor a bifurcation of fixed point in the correspond- 
ing classical system can give an quantitative account for the 
ground state, the fact that a superposition of two wave packet 
is a good approximation for the ground state is consistent with 
the classical observations of bifurcation. 

However, as shown in Figs.|2](c) and (f), the degree of en- 
tanglement between the spin and oscillator exhibits relatively 
larger deviations, especially for large SI. 

The poor accuracy of the variational solution in Eq. (|9]l 
is even more noticeable (more than 30%) on the momentum 
variance AP of the ground state. For both SI = 10 and 10^ 
case, the squeezing effect is considerably underestimated in 
this scheme. Moreover, for the true adiabatic regime, the vari- 
ational solution predicts that the squeezing effect suddenly 
disappears when g becomes larger than 1 . 

We can understand the reason for the failure by looking at 
the optimized variational parameters shown in Fig.|3](e). If 
g <C 1, the variational states are optimized at ai — —a2 ^ 
g^/ri/2. In this case, there is a considerable overlap between 
two coherent states which results in the squeezing effect. On 
the contrary, when g becomes larger than 1, the variational 
parameter abruptly changes, and the value for ai and a2 be- 
comes an order of 10 or 100 still having opposite signs with 
each other. Then, the overlap between two coherent states 
vanishes, and cannot give the squeezing effect. This observa- 
tion lead us to conclude that a deformation of individual wave 



packets is considerable in the true ground state, and unless 
the effect is taken into account in the trial wave function, the 
squeezing effect is inevitably underestimated. 

Motivated from the above observation, we propose a su- 
perposition of two displaced-squeezed states as a trial wave 
function, 

\'ip{ri,ai,r2,a2,t)) = (1 - i)|ri,ai) + i|r2,a2), (10) 

where \r, a) denotes the displaced-squeezed state ESl defined 
by |r, a) ^ S{r) \a) with S{r) = exp[(r6t6t _ r*bb)/2]. We 
refer to the wave function in Eq. ( [T0| as the double squeezed 
state (DSS). Without a loss of generality, the displacement 
parameter a and the squeezing parameter r are assumed to be 
real. 

The variational parameters ri,ai,r2,a2, and t are de- 
termined by minimizing the energy E{ri,ai,r2,a2,t) = 
{ip\H^ \tp) / (tplip) . It includes two direct terms. 



{rj,aj\H^\rj,aj) = sinhVj + (A - a^e'^f - ^e" 
for j = 1,2, and a cross term 
{ri,ai\Hi\r2,a2) 

— ; ■ r — — secn(ri — r2) smnri smnr2 

{ri,ai\r2,a2) 

+ [A — sech(ri — r2){ai coshr2 + a2 sinhri)] 
X [A — sech(ri — r2)(ai sinhr2 + a2 coshri)] 
1 



(11) 



where 



{ri,ai\r2,a2) 



-ilexp [—2 scch(ri — r2)aia2] (12) 



exp [aia2 sech(ri — r2)] 



^cosh(ri - r2) 



X exp 



(a? 



^l) + ("i - "2) tanh(ri - r2) 
2 



(13) 



The optimum variational parameters ri, r2, ai, Q!2, and t thus 
determined are plotted as a function of g in Fig. [3] The squeez- 
ing parameters ri and r2 of the individual wave packets in the 
DSS is indeed large in the regime g ~ 1, where the DCS vari- 
ational solution changes abruptly. This is consistent with our 
prediction. 

Using this variational solution, we calculate the momentum 
variance. The variance AP of the momentum P — i{W — h) is 
given by AP = (^/ilP^j-!/;) / (V'lV') in our case {ip\P\^) = 
0. The direct terms and cross term are respectively given by 



{rj,aj\P^\rj,a-j) 

(ri,Q;i|p2|r2,Q;2) 
(ri,ai|r2,a2) 



(j = l,2) and 



4(Q!ie''i -a2e''2)^ 



=2ri 



e 



2r2^2 



(14) 



Thus calculated momentum variation is shown in Fig. [2] (d) 
and (h). The squeezing effect is accurately captured by the 



DSS ( 10 1, and the accuracy can be attributed to the deforma- 
tion of the constituent wave packets. Our variational solution 
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FIG. 3: (color online) Optimized variational parameters for the dou- 
ble squeezed state [Eq. lj9}] for f2 = 10 (a,b,c) and for = 10^ 
(d,e,f). The variational parameters for the double coherent state 
[Eq. l |10^ ] are qualitatively similar (not show) except for ri = r2 = 
0. (a,d) The relative weight t in the superposition. (b,e) The dis- 
placement parameters a\ and 012. (c,f) The squeezing parameters r\ 
and r2- 

illustrates well the phase transition-like behavior for O — > 00 
regime. 

Interestingly, the optimal variational parameters satisfy 
ai sa — a2 and ri w r2 for ^ 1 [it is not the case for 



VI ^ 1; compare Fig. |3](b,c) with Fig. [3](e,f)]. In this case, 
the DSS in Eq. ( fTO] ), which has been expressed in terms of the 
spin-oscillator hybrid mode b, takes a simple form (not nor- 
malized) 

\ip{r, a, r, -a, t)) = (|r, a) + \r, -a)) <g) \i) 

+ (l-2<)(|r,a)-|r,-a))®|t) (15) 

in terms of the original spin and oscillator mode. This form 
shows transparently the entanglement feature in the optimal 
wave function: For 5<C 1, t« 1/2 and the spin and oscillator 
is separable. As g increases, t decreases, which leads to strong 
entanglement between spin and oscillator mode as shown in 
the insets of Fig.|2](c) and (g). and the overlap between \a) 
and |— a) {a ^ A) gives a finite amount of squeezing. 



III. DISCUSSION 

We have studied the spin-boson model ([T]i in the ultra- 
strong coupling limit by means of an effective classical po- 
tential and an improved variational wave function. 

We note that in the true adiabatic regime (il ^ 1), the crit- 
ical strength g for the entanglement coincide with the point 
where the squeezing effect is maximum. This coincidence 
can also be noticed from the crossover behavior of the effec- 
tive classical potential (|9]l. Nevertheless, the precise relation 
among the bifurcation in the classical potential in Fig. ([T]i or in 
the corresponding classical model |12], the squeezing effect, 
and the entanglement remains an interesting open problem in 
the ultra-strong coupling limit of the model ([T]i. 
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